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ABSTRACT 

Using WMAP 5 year data, we look for the average Sunyaev-Zel'dovich effect (SZE) 
signal from clusters of galaxies by stacking the regions around hundreds of known X- 
ray clusters. We detect the average SZE at a very high significance level. The average 
cluster signal is spatially resolved in the W band. This mean signal is compared with 
the expected signal from the same clusters calculated on the basis of archival ROSAT 
data. From the comparison we conclude that the observed SZE seems to be less than 
the expected signal derived from X-ray measurements when a standard /3-model is 
assumed for the gas distribution. This conclusion is model dependent. Our predictions 
depend mostly on the assumptions made about the core radius of clusters and the slope 
of the gas density profile. Models with steeper profiles are able to simultaneously fit 
both X-ray and WMAP data better than a /3-model. However, the agreement is not 
perfect and we find that it is still difficult to make the X-ray and SZE results agree. A 
model assuming point source contamination in SZE clusters renders a better fit to the 
one-dimensional SZE profiles thus suggesting that contamination from point sources 
could be contributing to a diminution of the SZE signal. Selecting a model that better 
fits both X-ray and WMAP data away from the very central region, we estimate the 
level of contamination and find that on average, the point source contamination is on 
the level of 16 mJy (at 41 GHz), 26 mJy (at 61 GHz) and 18 mJy (at 94 GHz). These 
estimated fluxes are marginally consistent with the estimated contamination derived 
from radio and infrared surveys thus suggesting that the combination of a steeper gas 
profile and the contribution from point sources allows us to consistently explain the 
X-ray emission and SZE in galaxy clusters as measured by ROSAT and WMAP. 
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1 INTRODUCTION 

In the last few years years, there has been an ongoing de- 
bate about the possibility that we do not yet quite under- 
stand the physics involved in the interaction between the 
gas inside galaxy clusters and the cosmic microwave back- 
ground (or CMB). This interaction between CMB photons 
and galaxy clusters is described by the Sunyaev-Zel'dovich 
effect (or SZE hereafter; Sunyaev & Zel'dovich 1972). The 
motivation for this debate is the apparent weaker strength of 
the SZE signal in galaxy clusters than expected from X-ray 
observations of the same clusters (Lieu et al. 2006, Afshordi 
et al. 2007, Atrio-Barandela et al. 2008). On the other hand, 
other works claim that there is no missing signal and that 
the interaction between galaxy clusters and the CMB can 
be well described in terms of the SZE once a realistic model 
for the gas distribution is considered (Mroczkowski et al. 
2009). Understanding the SZE is important also for cosmo- 
logical studies based on clusters and the SZE (see for in- 
stance Bartlett & Silk 1993, Barbosa et al. 1996, Aghanim 



et al. 1997, Diego et al 2002b, Carlstrom et al. 2002, Ko- 
matsu & Seljak 2002, Huffenberger et al. 2004). 

In the first attempt to detect the SZE in WMAP data, 
Diego et al (2003) found no SZE when correlating WMAP 
one year maps and maps derived from the ROSAT All Sky 
Survey (or RASS hereafter, Snowden et al. 1997). In this 
work the authors were not selecting only areas of the sky 
with known clusters; rather a whole section of the sky was 
used in the analysis. In more recent years, Atrio-Barandela 
et al. (2008) used a large sample of know galaxy clusters 
and observed that while the observed signal at small radii 
could be consistent with the expected SZE signal from the 
central parts of the clusters, the signal at larger radii (1 de- 
gree or more) is not consistent with the predictions. In the 
same work it was suggested that the pressure profile follows 
a Navarro-Frenk- White (or NEW, Navarro et al. 1996) pro- 
file rather than a /3-model profile. This result agrees with 
higher resolution X-ray observations where profiles steeper 
than a standard /3-model seem more adequate to fit the ob- 
servations (Neumann 2005, Vikhlinin et al. 2005, 2006, Cros- 
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ton 2008) and also with theoretical predictions (Hallman et 
al. 2007). In a diflterent work, Myers et al. (2004) used the 
APM Galaxy Survey and the Abell-Corwin-Ollowin (AGO) 
catalogue and through a correlation analysis with WMAP 
data, found evidence that significant amounts of gas lay in 
the outskirts of galaxy clusters (at distances up to 5 Mpc). 
Hernandez-Monteagudo & Rubino-Martin (2004) used an X- 
ray catalog of clusters to cross-correlate with WMAP first 
year data. They found a 2-5 a correlation signal and a mean 
temperature decrement of 15-35 )iK in the Raylcigh- Jeans 
region. However, no comparison was made with the expected 
SZE signal from these clusters. One of the most closely re- 
lated works to the present paper can be found in Lieu et al 
(2006) where they use a sample of 31 X-ray clusters and look 
for the SZE in WMAP one year data. They found evidence 
that there is less SZE than expected. The strength of their 
results is based on the fact that the models they used to pre- 
dict the SZE are consistent with the X-ray emission observed 
for each one of the 31 clusters. Afshordi et al. (2007) also 
found evidence of missing SZE signal using a sample of 193 
clusters. Our work will have some similarities with Lieu et 
al. (2006) but also some important differences, for instance 
the much larger size of our cluster sample. Also, we consider 
different models which could explain some of the differences 
between the results. The reader will also find some similari- 
ties with the work presented in Atrio-Barandela et al. (2008), 
in particular the cluster sample used in both papers is very 
similar but the analysis (and conclusions) differ significantly. 

In this paper we follow previous work and revisit the 
debate, but also perform a different analysis using the most 
accurate CMB data to date, namely the five year WMAP 
data (Hinsaw et al. 2009). We use a large sample of known 
X-ray clusters and stack the GMB data in an area around 
each one of them. The stacking procedure effectively reduces 
the contributions of instrumental noise and fluctuations in 
the CMB itself relative to the SZE making it possible to 
see the clusters (or their average signal) in WMAP data. 
For this process, it is crucial to have a large sample of clus- 
ters because otherwise the noise, and especially the GMB 
fluctuations, might not average to zero. Using the known 
positions of galaxy clusters offers the unique advantage of 
knowing were to look for the SZE. But it also provides us 
with a way of modeling the expected signal. Using the X-ray 
fluxes we are able to predict the SZE for each cluster (and 
hence, the average signal when the SZE signal of all clus- 
ters is added together) and compare the average predicted 
SZE with the observations. We find that the expected SZE 
signal depends strongly on the assumptions made about the 
internal distribution of the cluster gas. Our results neither 
strongly contradict nor strongly confirm the findings in some 
previous works where a lack of signal from clusters at mm 
wavelengths was found. As we will show, our results still 
show a lack of SZE signal in galaxy clusters but the dis- 
agreement between the observed signal and the expected 
one is significantly less than previously reported by Lieu et 
al (2006). 

We also explore some of the possible mechanisms re- 
sponsible for the apparently missing SZE signal. One of the 
most likely explanations for the lack of SZE signal is (apart 
from the internal distribution of the gas) contamination from 
radio and/or infrared galaxies inside the cluster. Since the 
SZE signal is negative at the frequencies used by WMAP and 



the radio and infrared emission is positive, a reduction of the 
SZE signal might occur from this contamination. Lin et al. 
(2009) made a survey of radio galaxies inside galaxy clusters 
in the frequency range 4.9-43 GHz. They found that half the 
observed sources had a step spectrum (as expected) but they 
also found that most of the unresolved sources had fiat or 
even inverted spectrum. These flat and inverted spectrum 
sources have more flux at higher frequencies than predicted 
from standard extrapolations (that usually assume a steep 
spectra) and they could be a significant source of contami- 
nation even at WMAP frequencies. On the other extreme of 
the frequency, SGUBA observations of galaxy clusters (Zem- 
cov et al. 2007) also reveal a significant amount of infrared 
emission in galaxy clusters at SCUBA frequencies. We ex- 
amine both in section 8. Wc discuss there whether the point 
sources are or are not significantly reducing the amount of 
SZE signal in WMAP. 



2 SUNYAEV-ZEL'DOVICH EFFECT AND 
X-RAYS IN CLUSTERS 

When CMB photons cross a galaxy cluster there is a prob- 
ability that some of the photons interact with the free elec- 
trons in the hot plasma through inverse Compton scatter- 
ing. Through this interaction, the GMB photons can gain 
a small amount of energy. This change of energy distorts 
the blackbody spectrum of the CMB photons. At frequen- 
cies smaller than 217 GHz the distortion results in a deficit 
of photons compared with the undistorted blackbody spec- 
trum. These missing photons appear at (or are promoted to) 
higher frequencies creating an excess in flux above 217 GHz. 
This is known as the Sunyaev-Zel'dovich effect (Sunyaev & 
Zel'dovich 1972). At WMAP frequencies, the clusters will 
appear colder than the CMB since we are in the frequency 
range where the deficit occurs {v < 217 GHz). The temper- 
ature decrement (increment for frequencies i> > 217 GHz) 
observed in the direction 9 can be described as; 

AT{e)=To j HeTdl (1) 

where To contains all the relevant constants including the 
frequency dependence {gx = x{e^ + l)/(e^ — 1) — 4 with 
X = hv/kT), Ue is the electron density and T is the elec- 
tron temperature. The integral is performed along the line 
of sight. An interesting quantity is the integrated tempera- 
ture decrement which can be obtained after integrating the 
above equation over the solid angle of the cluster (dfi) 

/M T 
n^Tdfldl oc 7° (2) 

where Mgas is the total gas mass of the cluster and Da{z) the 
angular diameter distance (instead of the gas mass one can 
use the total cluster mass if Mgas = ft * M and the baryon 
fraction, /(,, is constant). The advantage of using this expres- 
sion is that the total integrated signal is independent of the 
internal distribution of the gas. Once the mass, temperature 
and redshift of the cluster are know, a robust estimate of the 
integrated ATrot can be derived. The temperature can be 
obtained from X-ray observations of galaxy clusters. The 
mass can also be derived from the same observations using 
scaling relations between the X-ray luminosity and the total 
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Figure 1. y3-model (thin lines) compared with ADOS model (thick 
lines). In both cases, the dotted line represents the electron den- 
sity (ne), the dashed line the electron pressure {rigT or SZE) 
and the dotted-dashed line the X-ray emissivity (n^T"-^). The 
/3-model corresponds to an isothermal model with T = 3 keV, 
(3 = 2/3, r-c = 150 kpc, and rio = 0.01 cm"-''. The ADOS model 
has no = 0.01 cm'^, To = 3 keV, a = 3rc = 450 kpc, t = 0.2, 
a = 0.15 and n = 4. The ADOS model is explained in the text. 



mass. The total SZE flux is determined from equation (2) 
and the appropriate frequency dependence conversion be- 
tween brightness and temperature. The above expressions 
describe the thermal SZE. If the gas in the cluster is moving 
with a significant bulk velocity with respect to the CMB, 
another effect, the kinetic SZE, occurs in addition to the 
thermal SZE. In this work we will not consider the kinetic 
effect since we will be working with the average signal of a 
large number of clusters, and in this case the average bulk 
velocity can be neglected. 

The same electrons that interact with the CMB photons 
emit energy in the form of X-rays through the bremstrahlung 
process. An interesting situation is when the plasma has 
very high temperatures (T > 10 keV). In this case, a non- 
negligible fraction of the electrons have relativistic velocities 
(and energies) and are able to promote some of the CMB 
photons into the X-ray band itself. This mechanism increases 
the amount of observed X-rays and reduces the amount of 
CMB photons. This interaction has been suggested in the 
past as a possible explanation for the reduced SZE predicted 
from X-ray observations (Lieu & Quenby 2006). However, 
few of the clusters we consider have T > 10 keV. Considering 
only the non-relativistic case, the observed X-ray flux in the 
same direction 6 from a cluster is given by; 



Sxifi) — So 



JnlT^'^dl 



(3) 



where Di [z] is the luminosity distance. The quantity So con- 
tains all the relevant constants and corrections (including 
the gaunt factor, band correction and k-correction). 

An obvious advantage of combining X-ray and SZE ob- 
servations is that it is in principle possible to break the de- 
generacy between the multiple density profiles due to the 
difference dependence of the two signals on rig. 

One of the most common assumptions made when mod- 
eling the gas distribution in galaxy clusters is to assume that 
Tie is well described by a /3-model. 



[1 + (r/rc)2]3/3/2 

The /3-model is well motivated by observations of galaxy 
clusters and theoretical predictions. The exponent /3 typi- 
cally ranges from values of 0.5 to 0.9 but for practical rea- 
sons, a value of /3 = 2/3 is usually assumed (this choice 
allows for an analytical solution of the integral of ne). Other 
more elaborate (and accurate) models can be used instead. 
As an alternative model we will consider the following de- 
scribed in Ascasibar & Diego (2008). We will refer to this 
model as ADOS from now on. Ascasibar & Diego (2008) have 
shown how this model is capable of fitting high resolution 
profiles of X-ray clusters obtained with CHANDRA data 
(Vikhlinin et al. 2006). The ADOS model describes the ra- 
dial structure of clusters in terms of five free parameters: 
the total mass M of the system (or the chaxacteristic tem- 
perature To), a characteristic radius a, the cooling radius in 
units of a, < Q < 1, the central temperature in units of 
To, < t < 1, and the asymptotic baryon fraction at large 
radii, / ~ 1. Ascasibar & Diego (200S), have shown that 
this model is able to reproduce the observed X-ray proper- 
ties of real clusters. The total (gas plus dark matter) density 
follows a Hernquist (1990) profile. 



p(r) = 



M 



(5) 



r/a{l + r/aY ' 
whereas the gas temperature varies according to 

Tir) = ^" ^ + 

1+r/a 1 + r/ac' 

From hydrostatic equilibrium the mass-temperature re- 
lation can be obtained : 

(n + 1 = 

/imp a 

and the gas density profile is 



(6) 
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where n is an effective polytropic index, k is the Boltzmann 
constant, mp denotes the proton mass, ji ~ 0.6 is the molec- 
ular weight of the gas, and no is the central gas density. As 
in ADOS, we set n = 4 in order to obtain a constant baryon 
fraction at large radii. From Ascasibar & Diego (2008) we de- 
rive average values for the t and a parameters of the model, 
in particular we take t = 0.2 and a = 0.15. For the pa- 
rameter a we use a = 3.3rc which makes the ADOS model 
resemble a /3-model with the same rc at intermediate radii. 

A comparison of the density profiles of the two models 
is shown in figure 1. In the same figure we also show the SZE 
and X-ray emissivity for the two models. In our calculations, 
we assume that the /3-model is isothermal (T = const.). The 
ADOS model is steeper than the /3-model. This has impor- 
tant implications since the gas is more concentrated in the 
center and therefore less gas is needed in order to explain 
a given X-ray fiux. Also, dumpiness inside the cluster can 
create more X-rays with the same amount of gas due to the 
same mechanism (since the X-ray emissivity varies as n^). 
For instance, one can assume that some amount of dumpi- 
ness exist around the galaxies belonging to the cluster. High 
resolution observations with CHANDRA have reveled de- 
partures from the smooth and continuous X-ray emissivity 
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Figure 2. Luminosity temperature relation. The diamonds cor- 
respond to the BCS and eBCS catalog while the dots correspond 
to the REFLEX catalog. 

assumed in the construction of figure f. An example are 
the so called cold fronts first observed by CHANDRA (see 
for instance Markevitch et al. 2000, Vikhlinin et al. 2001, 
Mazzotta et al. 2001 and others). In the X-ray band, the 
cold fronts appear as a shock wave with significantly more 
surface brightness inside the shock wave than outside. The 
interpretation of these shock waves is that they are pro- 
duced by the oscillations of the gaseous core around the 
peak of gravitational potential. When the gas oscillates and 
moves to a region with lower potential energy, it expands 
and cools. Then, because it is cooler, it collapses to higher 
density, thus increasing ni and the X-ray luminosity (Mathis 
et al. 2002, Nagai & Kravtsov 2003, Ascasibar & Markevitch 
2006). Studies have also revealed that the pressure remains 
constant across this shock wave (e.g. Mazzotta et al. 2001, 
Nagai & Kravtsov 2003). This means that while the X-ray 
brightness can be enhanced by mechanisms like the cold 
fronts, the SZE remains constant, since the pressure does 
not change. We return to this issue in section 10 



3 X-RAY CLUSTER CATALOGS 

There are a variety of cluster catalogs. For the purpose of 
this work the most interesting ones are those based on an 
X-ray selection. There are two important advantages when 
using a cluster catalog to look for the SZE. The first one 
is obvious: knowing the positions of galaxy clusters allows 
us to select those portions of the sky where a cluster exists. 
Second, as mentioned earlier. X-rays are emitted by the same 
hot gas that produces the SZE. The X-ray signal of a galaxy 
cluster then allows to constrain the model used to describe 
the internal gas distribution. 

For this work, wo build a largo cluster catalog combining 
three catalogs: the Bright Cluster Sample (or BCS Ebeling 
et al. 1998), the extended Bright Cluster Sample (or eBCS 
Ebehng et al. 2000) and the ROSAT-ESO Flux-Limited X- 
Ray catalog (or REFLEX; Bohringer et al. 2004). These 
catalogs are X-ray selected using data from the ROSAT All 
Sky Survey (or RASS). The combined catalog covers both 
the north and south Galactic hemispheres for Galactic lati- 
tudes roughly above 20° and below -20° . There are 750 clus- 



ters in our combined catalog outside the Galactic exclusion 
zones. The combined catalog contains measured redshifts. 
X-ray fluxes, temperatures (for some clusters) and luminosi- 
ties for the brightest X-ray clusters observed with ROSAT. 
In the case of REFLEX clusters, the catalog docs not con- 
tain temperature estimates, but these can be derived from 
scaling relations as described below. A scaling relation was 
also used by Ebeling et al. (2000) to derive temperatures for 
some of the clusters in the eBCS catalog (see figure 2). 

Any model that we use to predict the SZE in clusters 

can be compared with the observed X-ray emission in those 
clusters. In an ideal scenario, we would like to do this ex- 
ercise on a one-by-one basis; Filing a model to the X-ray 
surface brightness of a particular cluster and using the same 
model to predict the amount of SZE in that cluster, then re- 
peating the process for every single cluster. This kind of ex- 
ercises has been done in the past (Lieu et al. 2006) but with 
a small number of very bright clusters. In section 10 wc will 
discuss in more detail the interesting findings of Lieu et al. 
(2006) . Only very bright clusters allow a reliable estimation 
of cluster profiles with ROSAT data. A small (but growing) 
number of clusters have been observed with CHANDRA and 
XMM (VikhUnin et al. 2006, Croston et al. 2008). These ob- 
servations allow a much more detailed analysis. 

In this work we adopt a different approach in which the 
sample size matters. We focus on having a large sample of 
clusters (even if their profiles arc not resolved) rather than 
a small but well defined sample of clusters. Instead of fitting 
the profile of a few dozen clusters we fit the average flux 
of a much wider sample of hundreds of clusters. Taking a 
large sample of clusters is crucial for reducing the CMB and 
instrumental noise fluctuations that contaminate the SZE 
signal. We use the X-ray observed flux in every cluster to 
normalize the gas model. By normalizing we mean select- 
ing the particular value of Uo in equations (4) and (8) such 
that the total flux predicted for our model coincides with 
the observed flux. Of course, this alone does not guarantee 
that one-by-one we are achieving a good fit to the surface 
brightness of every cluster. We are more interested in the 
average properties of clusters. For instance, most clusters 
present significant departures from circular symmetry but 
the average profile of a few clusters already starts to show 
this symmetry. 

Average properties can usually be described by scaling 
relations, and are powerful tools that have allowed us to set 
constraints on cosmological models using galaxy cluster ob- 
servations (Oukbir & Blanchard 1992, Lupino & Gioia 1995, 
Donahue 1996, Eke et al. 1996, Kitayama & Suto 1997, Ouk- 
bir & Blanchard 1997, Mathiesen & Evrard 1998, Donahue 
& Voit 1999, Diego et al. 2001, Henry 2004). Some of these 
scaling relations show very tight correlations between cluster 
properties. The best studied so far use X-ray data. In par- 
ticular, one of the most famous is the correlation between 
the X-ray luminosity [L^) and the plasma temperature (T) 
or Lx — T relation (David et al. 1993, Mushotzky & Scharf 
1997, Vikhlinin et al. 2002, Ettori et al. 2004). N-body sim- 
ulations, as well as analytical models, predict similar corre- 
lations between quantities derived from SZE observations. 
Scaling relations are representative of the average proper- 
ties of clusters and we use them to model the average sig- 
nal observed by two very different experiments, ROSAT and 
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Figure 3. Relation between tiie X-ray luminosities and the core 
radius. The diamonds represent the BCS and eBCS VTP radius 
(see text) divided by a factor 3 and the dots represent the rough 
estimate of the core radius given in the REFLEX catalog. The 
solid line shows the theoretical relation expected for a self-similar 
model while the dashed line represents and alternative and arbi- 
trary model with slope 0.4. 



WMAP. Two examples of scaling relations used in this work 
are shown iu figures 2 and 3. 

In figure 2 we show the Lx — T relation for clusters in 
the different catalogs. The diamonds show the corresponding 
temperature and X-ray luminosity (in the 0.1-2.4 keV band) 
for BCS and eBCS clusters while the dots correspond to the 
REFLEX clusters. In the case of BCS and eBCS clusters, 
the X-ray luminosities were given originally for a different 
cosmological model than the one used in the REFLEX cat- 
alog. Starting from the fluxes, we compute the luminosities 
for BCS and eBCS clusters but for the same cosmological 
model used in the REFLEX catalog. Some of the clusters in 
the BCS and eBCS catalogs did not have temperature es- 
timates. In these cases, the temperature was obtained from 
the luminosity using a scaling relation. In the case of RE- 
FLEX clusters, none of the clusters in the catalog had tem- 
perature estimates. We used a scaling relation similar to the 
one in the BCS and eBCS clusters. 



Lx = 1. 



erg/s, 



(9) 



where L^ is the luminosity in the 0.1 — 2.4 keV band and T 
is the gas temperature in keV. 

A relation between the X-ray luminosity and physical 
core size of the cluster can be obtained from the above Lx~T 
relation and the following two scaling relations for a model 
that we will refer as self-similax model. These are the rela- 
tion between the virialized mass and temperature (Horner 
& Mushotzky 1999, Diego et al. 2001) 



! 8M 



2/3 



keV 



(10) 



and a relation between the virialized mass and the virial 
radius. (Horner & Mushotzky 1999, Diego et al. 2001) 



r„ ^ 1.3M^/-^(1 + 2)"^ Mpc 



(11) 



In the above expressions, M is the mass in units of 10^^ Mq. 
If we assume a fixed ratio of the virial to core radius r„/rc = 
p then we finally obtain. 
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kpc 



(12) 



where Lx is the luminosity in the 0.1 — 2.4 keV band in units 
of 10*'* erg/s. This scaling law is shown in figure 3 (solid 
line) and compared with a rough estimate of the core radius 
for REFLEX clusters (dots) as given by Bohringer et al. 
(2004). In the case of the BCS and eBCS clusters, we plot the 
Voronoi Tesellation Percolation (VTP) radius, divided by 
an arbitrary factor of three to omake the BCS-|-eBCS radii 
consistent with the REFLEX estimates. The VTP radius 
does not have a physical meaning but rather is the radius 
at which the cluster is seen as an overdense region in the 
ROSAT PSPC detector with the VTP algorithm (isophotal 
radius). The diamonds in figure 3 should not be confused 
with the core radius but rather with a very rough estimate 
of the physical size of the core region in the BCS and eBCS 
clusters. Scaling relations similar to the one shown in figure 
3 have been found in the past between the isophotal size and 
the X-ray temperature (Mohr et al. 2000). Also shown in the 
same figure is a different (and arbitrary) model where the 
core radius scales as I/JJ '*. We will use this arbitrary model 
as a way to test significant departures from the self-similar 
model. Observationally, the radius-luminosity relation is not 
as well established as the Lx—T relation. As we will see later, 
this will be probably the most important limiting factor on 
our conclusions as the predicted signal depends strongly on 
the core radius. In particular, in both cases. X-rays and SZE, 
the total integrated flux scales as r^. This means that a 
factor 2 difference in the scaling relations leads to an order 
magnitude difference in the predicted flux. The good news is 
that the combination of X-ray and SZE measurements helps 
reduce the range of possible models. 

An interesting feature of our normalization process is 
that since we require that the X-ray flux from our model 
match the observed flux, smaller core radii will boost the 
central electron density (which enters like rio in the predicted 
X-ray flux) in order to compensate for the rl decrease. As 
a consequence the predicted cluster signal will appear more 
concentrated in the center but with a large amplitude rela- 
tive to models with larger core radii. Since the X-ray tails of 
the gas density are difficult to see (or detected by an X-ray 
experiment) the cluster will look as if it is brighter when 
taking smaller core radii (this is of course a visual artifact 
as the total flux is the same for all models given a particular 
cluster). In contrast, the predicted flux in the SZE case is 
not constrained, and can be larger or smaller depending on 
the model. It is instructive to look at the order zero expres- 
sion for the expected total fluxes. In the case of the SZE, the 
integrated total flux is proportional to Uo x (times a slow 
varying function which depends on the ratio p = fy/rc). On 
the other hand, the total X-ray flux goes like x (again 
times another geometric factor with a weak dependence on 
p). Since the total X-ray flux is constrained to a value. Sobs, 
the electron density is just Uo = x rj^^^. Substitution 
of Uo into the SZE flux results in the SZE having a total 



flux going like S^^ 



1/2 



3/2 



for a given cluster. That is, while 



the total X-ray flux is constrained to be equal to Soba, the 
total predicted SZE flux is not constrained and grows with 
the core radius as r?''^. 
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Figure 4. Average cluster profile as derived from RASS maps. 

The fields of view in this figure and figure 6 is 2.51° X 2.51° and 
both arc centered on the clusters in the X-ray catalog. Also shown 
(on the same scale) are the PSF of the ROSAT PSPC detector 
(bottom left), the PSF introduced by the HEALPIX pixelization 
(bottom middle) and the PSF due to RASS 12 arcminute pix- 
elization (bottom right). 



4 STACKING OF ASTROPHYSICAL DATA 

When looking at individual clusters in the WMAP data, 
their signal to noise ratio is too small to allow a clear de- 
tection of most individual clusters. The noise in this case 
is a combination of the instrumental noise (well approxi- 
mated to first order as a random Gaussian noise) plus the 
CMB. At a lower level. Galactic foregrounds and possibly 
unresolved extragalactic sources (i.e radio and/or infrared 
galaxies) also contribute to the noise fluctuations. One way 
to increase the signal to noise ratio for galaxy clusters is by 
stacking the area around known galaxy clusters (observed 
for instance by their X-ray emission). The stacking proce- 
dure reduces the contribution of the CMB and instrumental 
noise since these two signals should average to zero for a 
sufficiently large number of stacked fields. In the case of the 
Galactic foregrounds and extragalactic unresolved sources, 
these two components do not average to zero in the stacking 
process, but will contribute as a largely uniform background 
to the average map (especially if one considers areas of the 
sky far enough from the Galactic plane). 

In our stacking procedure we first construct full sky 
maps of the data to be stacked (in those cases where the 
all sky maps are not already available). An area (or field of 
view, FOV hereafter) of 2.5° x 2.5° is then selected around 
each galaxy cluster in our combined catalog and the area 
is projected into a plane using an orthographic projection. 
We sample the FOV of 2.5 degrees with 256 x 256 pixels 
(i.e a pixel size of « 0.6 arcmin). Due to the random ori- 
entations introduced by the stacking procedure, the original 
pixels in the all-sky maps take random orientations and do 
not align with other pixels in different FOVs. The stacking 
procedure introduces an average pixel resulting from the av- 
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Figure 5. One dimensional profile of RASS stacked clusters (solid 
line) compared with the normalized effective point spread func- 
tion (or psf, dotted line). This psf contains the combined effect of 
the intrinsic psf of the ROSAT pspc detector, the 12 arcminute 
RASS pixelization and the Nside = 1024 HEALPIX pixelization. 

eragc of all the possible orientations of the original pixel. In 
the next section we will discuss how to estimate this effect 
quantitatively, using the X-ray data as an example. 



5 X-RAY DATA/ROSAT DATA 

We use the diffuse X-ray maps derived from the ROSAT All 

Sky Survey (RASS) maps. These maps are described in de- 
tail in Snowden et al. (1997). The difi'use X-ray RASS maps 
are given in six Zenithal Equal Area (or ZEA) projections 
and with a pixel size of 12 arcminutes. 

The brightest point sources have been removed using 
the full resolution images of the PSPC detector. The ex- 
traction radius for these sources is at least 3.5 arcminutes. 
As a consequence of the bright source removal, the central 
peak in some bright clusters is also removed. 

The bright source catalog used in Snowden et al. (1997) 
to mask out point sources included at least 400 cluster can- 
didates (see a more detailed description of the bright source 
catalog in Voges et al. 1999). Some of those candidates 
correspond to real clusters and some not. We have cross- 
correlated our catalog containing 750 clusters and the bright 
source catalog of Voges et al (1999) and found that 314 of 
the objects in the Voges et al. (1999) bright source catalog 
satisfy the criteria of being at an angular distance of less 
than 1 arcminute from one of the clusters in our sample and 
have an extent of 30 arcseconds or more. Based on this es- 
timation we can say that we expect 40% of our sample to 
suffer from masking of their central region. We return to this 
issue in section 7. 

The RASS maps are given in different energy bands. We 
use the R6 band (0.91-1.31 keV) for which cluster emission 
is supposed to be maximal compared with emission from the 
Galaxy (dominant at lower energies) and the unresolved X- 
ray background (more relevant at higher energies). At each 
band the maps are presented in 6 azimuthal projections. 
We combine these 6 projections into a single all-sky map 
using the HEALPIX pixelization (Gorski et al. 2005). The 
original 12 arcminute RASS pixels are oversampled using 
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a 3.4 arcminute HEALPIX pixel (Nside=1024 in HEALPIX). 
We also subdivide each RASS pixel into 121 square pixels to 
avoid missing some of the smaller HEALPIX pixels. The final 
HEALPIX map contains over 12 million pixels and with an 
effective resolution which is the combination of 3 smoothing 
kernels. The first smoothing kernel represents the intrinsic 
Point Spread Function (or PSF hereafter) of ROSAT's PSPC 
detector. This PSF depends on the energy and the off-axis 
distance (Hasinger et al. 1994). We assumed an energy of 1 
keV (consistent with the energy range of the R6 band) and 
an off-axis distance of 20 arcminutes which is approximately 
the average distance of a serendipitous source in the field of 
view of the PSPC detector. With these parameters the PSF 
of ROSAT has a typical scale of 1-2 arcminutes (Hasinger 
et al. 1994). The second kernel represents the stacking pro- 
cedure of the original 12 arcminute RASS pixels. Since the 
center of a galaxy cluster can fall in any position within 
the 12 arcmirmtc RASS pixel (and not just in the center), 
this has to be taken into account. The stacking procedure 
averages all the possible orientations of the 12 arcminute 
pixel containing the center of each galaxy cluster. The spe- 
cific form of this kernel is given below. The third kernel is 
similar to the second one but it accounts for the stacking 
of the 3.8 arcminute HEALPIX pixels. Both, the second and 
the third smoothing kernel can be described by the following 
formula 

max 

) if < Omax (0 otherwise) (13) 

where is the angular separation to the center and Omax = 
l^Opix * \/2 with l^Opix being the pixel size (12 arcminutes 
for ROSAT and 3.7 arcminutes for HEALPIX) .This expres- 
sion results from fitting random orientations of pixels of size 
A.Opix and containing the central position of the source. This 
position can be anywhere (i.e is also a random variable) 
within the pixel. 

The final average map obtained after stacking the ar- 
eas 2.5° X 2.5° around the 750 clusters is shown in figure 
4. The angular extent of the average cluster is very simi- 
lar to the kernel representing the RASS pixel plus stacking 
process which is not surprising given the big pixel size of 
the RASS pixel (12 arcminutes) making most of the clus- 
ters appear as unresolved or barely resolved. Figure 5 shows 
the one-dimensional profile of the stacked X-ray fields and a 
comparison with the effective PSF resulting from combining 
the intrinsic PSF of ROSAT and the two smoothing kernels 
due to the stacking-l-pixelization process. 



6 WMAP DATA 

We use three of the WMAP frequency bands, Q (at 41 GHz), 
V (at 61 GHz) and W (at 94 GHz) to search for the signal 
from known galaxy clusters in the average stacked map. A 
more detailed description of the WMAP five year data and 
its main results can be found in Hinsaw et al. (2009) . Due to 
beam dilution effects, we expect to have the strongest SZE 
signal in the W band which has superior angular resolution. 
The average beam profile of WMAP is well approximated 
by a Gaussian in all three bands with FWHM — 0.51° for 
Q, FWHM = 0.35° for V, and FWHM = 0.22° for W. In 
addition to this beam, the pixelization-l-stacking procedure 
introduces an additional smearing function as described in 




Figure 6. SZ in WMAP data. Prom top to bottom, W, V and Q 
bands. Each panel shows the average signal obtained after stack- 
ing an area of 2.51° x 2.51° around the clusters in the X-ray cat- 
alog. Those clusters with bright radio sources have been masked. 
Note the very clear decrement in the W band. The decrement can 
still be seeing in the V and Q bands but with less intensity due 
to the beam smearing. 

section 5. The combination of the intrinsic WMAP beam and 
the pbcelization-hstacking procedure produces a slightly dif- 
ferent effective beam (see below) . In order to remove bright 
radio sources, wc use the kpO mask. 

Following a process similar to the one described in the 
previous sections, we stack the 750 fields around the X-ray 
clusters and for each band {Q, V and W). The result is 
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Figure 7. One dimensional average cluster profiles (dotted lines) compared with tlie effective beam profile (dashed lines). The effective 
beam profile is a convolution of the intrinsic WMAP beam and the pixelization+stacking smearing function. The effective beam has been 
normalized to the maximum amplitude of the average cluster profile. The x and y axis grow as the square of the scale and temperature 
respectively to better show the details in the central part and tails. For instance, the maximum decrement is ax;hieved in the W band 
with a amplitude of » 4.5^ = 20 ^xK (on the same plot 16 //K = 4^). Note how the average cluster profile seems to be resolved in the W 
and V bands. 



shown in figure 6. A clear signal is seen in the three bands 
at the position of the clusters. As expected, the signal is 
stronger in the W band, mostly due to its better angular 
resolution. Around the central position of the clusters we 
can see a degree scale structure in the three bands with hot- 
ter and colder fiuctuations. This is most likely the result of 
stacking the CMB arusotropies in the 750 fields of view. We 
tested this hypothesis by stacking the same immber of fields 
of view of a simulated map of the CMB. This simulation 
was done with the same angular resolution as the WMAP 
W band. The resulting map shows similar patterns to the 
ones observed in figure 6 and with the same strength. The 
average of the stacked simulated CMB was 4.8 yuK and the 
stacked map had a dispersion of 2.3 pK, similar to the ob- 
served dispersion in the borders of the maps in figure 6. As 
an additional test we looked at the stacking of the WMAP 
fields with the 80 largest (most positive fields or MPF) and 
80 smallest (most negative fields or MNF) values in a central 
region of 7 arcminute radius. The motivation for this exer- 
cise is to look for possible point source contamination in the 
stacked MPF. The same exercise was made with the simu- 
lation of the CMB described above (with no instrumental 
noise). In the case of the CMB simulation, the stacked MPF 
shows a large scale positive fiuctuation with a scale of 1-2 de- 
grees and an amplitude of about 100 ^K. The stacked MNF 
shows a similar structure but with opposite sign. When look- 
ing at the WMAP data, the stacked MPF shows a structure 
similar to the one seen with the simulated CMB. No signs 
of point source contamination are observed in these fields 
indicating that if point source contamination occurs it must 
be sub-dominant with respect to the CMB. In the case of 
the stacked MNF a similar structure but with negative sign 
was observed. However, in this case a small scale negative 
fluctuation is seen at the center of the field of view, specially 



in the W band indicating that among the sub-sample of 80 
MNF, there is a significant contribution from the SZE. 

When looking at the one dimensional profiles (obtained 
after averaging in rings around the center), we cam appre- 
ciate in figure 7 how a comparison with the effective beam 
indicates that the average cluster profile is resolved in the W 
band as well as in the V band. As noted above, the dispersion 
in the wings of the one dimensional profile is consistent with 
the 2.3 /iK dispersion of the simulated (and stacked) CMB 
maps. The profile observed in the lower resolution Q band is 
compatible with the profile (dashed line) expected for an un- 
resolved source convolved with the combined PSF. Resolving 
the average cluster profile with WMAP data means that ei- 
ther the average signal is dominated by low redshift clusters 
or that there is significant gas in the outskirts of galaxy clus- 
ters (a third option is of course a combination of the two). 
The second option is interesting since if confirmed it could 
be an indication that a fraction of the missing baryons could 
be in these regions. 

Another interesting result is the fact that if we do not 
use the kpO mask, and hence leave the bright microwave 
sources in the WMAP fields, the SZE signal disappears in 
the Q band and is significantly reduced in the V band while 
remaining basically unchanged in the W band. This sug- 
gests that significant radio source contamination from non- 
subtracted radio sources might still occur in the lower fre- 
quency bands while it should be much smaller in the W 
band. This, however, does not rule other kinds of point 
source contamination (like infrared sources), as the point 
sources on which the kpO mask is based are mostly (if not 
all) radio sources. 

We should note that the average of the stacked fields is 
not zero. Because we are using the uncleaned maps, there 
is still a non-negligible and positive contribution from the 
Galaxy. Before comparing the SZE profiles with models 
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Figure 8. Predicted versus observed (solid line) one dimensional 
profile for RASS clusters. The thick lines represent models that 
give a relatively good fit to the WMAP data while the thin lines 
are used for models who give a better fit to RASS data. The 
models represented as thin lines correspond to (dotted line) the 
standard /3- model (with 3 = 2/3), concentration parameter p = 
Tv/tc = ^ and a core radius equal to two times the one predicted 
by the self similar model. m This is denoted model Mi in table 
1. The thin dashed line is for M2 in the same table. The thin 
dot-dashed line is for M3. The thicker lines correspond to (dotted 
line) M4, thick dashed line M5 and thick dot-dashed line Mq. 

is necessary to subtract the Galactic contribution to the 
stacked fields. Since this value is unknown, wc subtract a 
constant from each map such that the one dimensional pro- 
files average to zero in the tails of the one dimensional pro- 
files (see figure 7). 



7 MODEL COMPARISON 

In this section we compare different models with the X-ray 
average emission from RASS maps and the SZE emission in 
WMAP data. It is important to note that in constructing 
the RASS average profile, wo have masked out the central 
part of the 314 clusters identified in the Voges et al (1999) 
catalog. The masking does not affect the computation of 
the central density, since the X-ray fiuxes were derived us- 
ing the full resolution PSPC data of Snowden et al. (1997) 
and therefore do not mask the central regions of the maps. 
However, we should not expect a perfect match between the 
predicted and observed average signal at the cluster cen- 
ter due to very non-linear processes such as cooling fiows 
that are not included in our model. In fact, due to these 
non-thermal effects, we should expect a good model to pre- 
dict slightly lower X-ray fiux than the observed one since 
we are assuming that all X-ray fiux is due solely to ther- 
mal brcmstrahlung (and emission lines). The wings of the 
average cluster profile should not be affected as much by 
these highly non-linear phenomena and hence they should 
be more useful to set reliable limits on the models. Any re- 
alistic model that we consider should predict correctly (or 
even better slightly undcrprcdict) the average signal seen in 
the RASS maps at radii between ~ 15 — 40 arcmin. 

Modeling the average X-ray emission from clusters as 
seen by ROS AT PSPC detector involves a series of steps that 
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Figure 9. Same as in figure 8 but showing the difference of the 
predicted minus observed one dimensional profiles for RASS clus- 
ters. 



Table 1. All /3- models have /3 = 2/3. SS stands for Self- Similar 
model and Alt for Alternative. The third column (rc/a) shows the 
core radius (for the /3-model) or the scale factor a (for the ADOS 
model). SSx4 means the self similar normalization has been multi- 
plied by a factor four and Altxl.5 means that the normalization of 
the Alternative model has been multiplied by a factor 1.5. SSxlO 
in the ADOS model means that the scale factor a is equal to ten 
times the self-similar core radius model. 



Model 


Type 


rc or a 


p 


Fits 


Ml 




SSxl.5 


5 


X-ray 


M2 


n 


Altxl.5 


10 


X-ray 


M3 


ADOS 


SSxlO 


5 


X-ray 


M4 


/3 


SS 


5 


WMAP 


Mb 




Alt 


5 


WMAP 


Me 


ADOS 


SSxia 


10 


WMAP 


Mr 


ADOS 


ssxia 


9 


Both ? 



need to be done carefully. First, we need a model for the gas 
distribution in throo-dimonsions (or 3D) that can be inte- 
grated to compute the X-ray luminosity and flux. The mod- 
els considered in this work consist of the standard /3-model 
and the polytropic model ADOS as described earlier. For the 
X-ray emission we use a Raymond-Smith model (Raymond 
& Smith 1977) with a metallicity of 0.3. This model includes 
emission lines which can change significantly the X-ray fiux 
for lower temperature clusters (T of few keV or less). The 
total X-ray flux is then converted into observed fiux in the 
0.1-2.4 keV band. This predicted flux can be compared with 
the observed one. This process involves two corrections, flrst 
the band correction that transforms bolometric flux into fiux 
in a given band (in the rest frame of the cluster) and second 
the i^— correction that accounts for the appropriate rcdshift 
distortions. For these corrections we used the values given 
in tables 3 and 5 in Bohringer et al. (2004) . In order to com- 
pare with the average profile obtained from RASS maps, 
the flux needs to be converted into counts per second (or 
cps) in the R6 band. For this conversion wo use the code 
XSPEC (Arnaud 1996) to fold the same Raymond-Smith 
model with the response matrix of the PSPC detector and 
in the appropriate range of energies of the RASS R6 band 
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(0.9-1.3 keV). The conversion factor depends on the plasma 
temperature and cluster redshift. For temperatures above 
a few keV and low redshift this conversion factor is about 
1 cps = 2 X lO"'^^ erg/s/cm^. We have assumed a negli- 
gible hydrogen column density (most clusters see less than 
1.5 X 10^^ cm~^ uh)- The error introduced by neglecting the 
hydrogen absorption is small (~ 10%) for the range of tem- 
peratures, energies and column densities used in this work. 

Given a particular assumption for the gas distribution 
in a cluster at redshift z, a model for the core radius rc 
and the central electron density rio, we can compute the ex- 
pected X-ray flux and SZE by integrating equations (1) and 
(3). The integration is carried out up to a maximum radius 
Rmax- The ratio between this maximum radius and the core 
radius is denoted by the parameter p and both the X-ray and 
SZE signals de4pend on the value of p although significantly 
less than on the choice for rc. The 3D models are projected 
along the line of sight and we obtain two-dimensional (or 
2D) maps of the SZE and the X-ray emission. The inte- 
grated X-ray flux is compared with the observed X-ray flux 
for each cluster. Since the integral in equation (3) is propor- 
tional to rio, our initial model can take any arbitrary value 
for no and then be re-scaled by the appropriate factor that 
makes the predicted X-ray flux match the observed one. The 
root-square of this scaling factor is the multiplied by the SZE 
emissivity to produce the predicted SZE signal for that clus- 
ter. The 2D maps are then stacked together and an average 
profile can be obtained for the SZE and X-rays. For build- 
ing the X-ray average profile we first convert the 2D fiux 
maps into their corresponding 2D cps maps in the 0.9-1.3 
keV band as described above. This conversion is done before 
stacking the maps since the conversion factor depends on the 
redshift and temperature of the clusters. As mentioned ear- 
lier, we remove the central 3.5 arcmin radius region in those 
clusters that were identified in the BSC catalog of Voges et 
al. (1999). After convolving the 2D maps with the appro- 
priate efi'ective PSF we can compare our predicted signal 
with the observed one in ROSAT and WMAP. The effective 
PSF for ROSAT was described in section 5. For WMAP 
the effective PSF is again the combination of the intrinsic 
WMAP beam and the smearing introduced by the pixeliza- 
tion-|-stacking procedure. This smearing function can be de- 
scribed again by equation (13) but with AOpix w 6.7 arcmin 
which is the pixel size used for WMAP data (see equation 
13 and description afterwards). 

In figure 8 we show different models that seem to agree 
well with the average X-ray profile. These models are listed 
in table 1. We have chosen two sets of models. The first set 
(represented with thin lines in figure 8) reproduce particu- 
larly well the average X-ray profile observed in RASS maps. 
The set with thicker lines still reproduce more or less this 
profile but seems to deviate in the wings, predicting less X- 
rays than observed. As we will see later, these models (thick 
lines) will fit the SZE WMAP profiles better than the thin 
line models. For a better comparison, the same models are 
represented again in figure 9 but this time showing the dif- 
ference between the observed profile and the predicted ones. 
When the same models are compared with the SZE we find 
that it is nearly impossible to find a model that can fit both 
data sets simultaneously (see figures 10 and 11). Models Mi, 
M2 and M3 give a good fit to the X-ray data but a bad fit to 
the SZE. Conversely, models M4, M5 and Me give a better 
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Figure 10. Models from figures 8 and 9 which fit the X-ray 
data well compared with WMAP data. Note how the /3-modeIs 
(dotted and dashed curves) over-predict significantly the observed 
WMAP average decrements. The models correspond to Mi (dot- 
ted), M2 (dashed), and M3 (dot-dashed) from table 1. The ver- 
tical axis shows the negative SZE decrement with opposite sign. 
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Figure 11. SZE fitting models from figures 8 and 9 compared 

with WMAP data. The models correspond to M4 (dotted), Mb 
(dashed), and Mg (dot-dashed) from table 1. 



fit to the SZE data but a poorer fit to the X-ray data. We ob- 
served that large core radii are needed to fit the X-ray data. 
As described earlier, the SZE grows as r?''^ when the gas 
profile is constrained with the X-ray fluxes. Consequently, 
the large Vc needed to fit the X-rays predict a large amount 
of SZE. Smaller values for Vc push the SZE down to levels 
compatible with the SZE seen in WMAP but this smaller rc 
would fail in fitting the wings of the X-ray profile. It seems 
that some level of compromise will be needed to reconcile 
both data sets. Given the assumptions made in the model- 
ing of the X-ray average signal, it is reasonable to assume 
that we do not expect a perfect match between the observed 
X-ray average profile and the predicted one in our model. 
For instance, no cooling flows or shock fronts were included 
in our model. In table 1 we list one model (A/7) that makes a 
relatively good fit to both data sets (see figure 16 below) but 
is not the best fitting model to either of the two data sets in- 
dividually. The model My under-predicts slightly the X-ray 
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emission over most of the radial profile. This, however, is a 
desirable feature since we should not expect a perfect fit for 
a model that does not include some phenomena like cooling 
flows, shock waves in the gas, or cold fronts. These mecha- 
nisms boost the X-ray emission without requiring more gas. 
Instead, a good model should slightly under-predict the X- 
ray emission in order to leave some room for these highly 
non-linear processes to fill the gap. The same processes are 
expected to have little impact on the SZE. In the next sec- 
tion we will sec that this model has another very interesting 
feature. It is reasonably consistent with estimates of the ex- 
pected level of contamination from point sources in galaxy 
clusters. This model can be considered as a concordance 
model in the sense that it delivers a reasonable fit consis- 
tent with the available observations. Note that no /3-model 
does as well. 



8 POSSIBLE CONTAMINATION FROM 
RADIO AND/OR INFRARED SOURCES 

Previous papers have reported that the amount of SZ ob- 
served in WMAP data is less than expected Atrio-Barandela 
et al. (2008), Lieu et al. (2006), Afshordi et al. (2007), Bielby 
& Shanks (2007). We confirm this result by looking at the 
temperature fluctuations in clusters. Also, the lack of signal 
seems to be more prominent in the outskirts of the clus- 
ter suggesting that the pressure proflle falls faster than pre- 
dicted by the /3- model (Atrio-Barandela et al. 2008) . A pos- 
sible mechanism to explain both observations would be a 
spatially uniform distribution of radio galaxies inside galaxy 
clusters. Such a uniform distribution would make the inte- 
grated radio signal more relevant at larger radii and also 
the overall positive contribution would reduce the total in- 
tegrated SZ signal. We check this hypothesis using a catalog 
of radio sources (sec below). It is known that some level 
of radio emission exist in galaxy clusters at the frequencies 
of WMAP (see discussion of the results of Lin et al. 2009 
below). In an earlier work, Cooray et al (1998) made a sur- 
vey of radio galaxies in clusters at 28.5 GHz. They found 
that there were 4 to 7 times more sources than predicted 
from a low-frequency survey in areas outside gala^xy clus- 
ters. They concluded that contamination from radio sources 
could be a serious threat to mm or cm SZE surveys. In a 
later work. Coble et al (2007) confirmed these results (also 
at 28.5 GHz) with a larger sample of sources in galaxy clus- 
ters. Lin et al. (2009) made an even more extensive study 
including also higher frequency bands. One of their most rel- 
evant conclusions for this work is that a significant fraction 
of radio sources show inverted or flat spectrum at frequen- 
cies beyond 10 GHz. The resuhs of Lin et al. (2009) will be 
discussed in more detail later. White & Majumdar (2004) 
predicted the average level of contamination of radio and 
infrared sources at mm wavelengths based on extrapolations 
from existing data. In particular, they predict that in the W 
band the confusion noise from infrared sources should equal 
the confusion noise from radio sources (with both equal to 
approximately 10 muK for a one arcminute beam). Other 
works predict the contribution of point sources based on an- 
alytical models (Toffolatti et al. 1998, Be Zotti et al. 2005). 
These works do not focus on galaxy clusters but rather on 
the average contribution of point sources to the CMB. Nev- 



ertheless they are still useful to set a minimum limit to the 
contribution of point sources to galaxy clusters (though we 
expect an enhanced signal due to clustering). 

At WMAP frequencies, we should expect the calcu- 
lated SZE signal to over-estimate the observed SZE by some 
amount, the difference being due to point source contamina- 
tion. In this section we try to estimate this level of contami- 
nation and combine it with the concordance model described 
above (M7). 

8.1 NVSS data 

The NRAO-VLA Sky Survey (or NVSS, Condon et al. 1998) 
catalog contains about 1.8 million sources with measured 
fluxes at 1.4 GHz. The NVSS survey has an effective reso- 
lution of 45 arcseconds, much smaller than the angular res- 
olutions relevant for this work and its effect can be safely 
neglected for our purposes. Using this catalog we build an 
all sky map using again the HEALPIX pixelization and the 
same pixel size as in the WMAP data (6.8 arcmin pixel). 
The stacking process at the positions of the 750 clusters on 
our list is similax to the cases discussed above. 

In figure 12 we show the result of the stacking proce- 
dure. NVSS bright sources can be seen at the center of the 
clusters in our catalog. We have divided the NVSS catalog 
into different flux bins and found that the NVSS sources 
seem to be strongly correlated with clusters down to fluxes 
of a few tens of mjy at 1.4 GHz. How much of this flux 
still remains at WMAP frequencies is still an open question 
but works like the ones presented in Lin et al. (2009) are 
helping to throw some light on this issue. From our results, 
however, we see no evidence of a spatially uniform distribu- 
tion of radio sources inside the clusters. Instead, the radio 
sources seem to be concentrated in the center of the clusters 
(as found by Lin et al. (2009)). 

8.2 Infrared data 

On the other side of the spectrum, we repeated the same 
exercise but on the infrared side using the IRAS galaxy 
catalog from Rowan-Robinson et al. (1991). This catalog 
contains 17711 sources with measured fluxes at 12 /um, 25 
/jm, 60 /im, and 100 jxra. For the stacking we used the 100 
pm fluxes as the corresponding frequency is closest to the 
WMAP ones. We found no signiflcant flux overdensity in 
the stacked position of our clusters. This is not surprising 
given the small redshifts {z < 0.02) of the galaxies in the 
IRAS catalog. Also, the large difference between the IRAS 
largest wavelength (100 fim) and the WMAP ones makes 
it difficult to extrapolate the fluxes down to WMAP fre- 
quencies. A related and interesting work can be found in 
Zemcov et al (2007). Using SCUBA archival data around 
galaxy clusters, the authors try to recover the SZE signal 
at 850 /im. Their sample contains 44 clusters. Among their 
conclusions, the most revealing are the fact that a large frac- 
tion of their sample shows significant emission from active 
galactic nuclei (or AGN) at the center of the clusters and 
that for the rest, the derived amplitudes of the SZE incre- 
ment are consistently higher than the amplitudes inferred 
from lower-frequency measurements of the SZ decrement. 
The disagreement could be easily explained if additional (i.e 
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Figure 12. Stacked NVSS sources at the positions of the clusters. Prom left to right, sources with fluxes smaller than 5 Jy, fluxes, smaller 
than 1 Jy, 0.5 Jy and 50 Mjy respectively. Also shown is the smearing introduced by the stacking+pixelization. 



not accounted for) emission from the infrared galaxies in the 
cluster is still contributing substantially to the total flux. In 
Zemcov et al. (2007), the authors provided also a very inter- 
esting list of 83 candidate sources detected in galaxy clusters 
together with their fluxes at 850 /itm. Fluxes in this sample 
range from 3 mJy to 70 mJy and with an average flux of 8.7 
mJy. Assuming a law for the flux (see for instance Greeve 
et al 2008) to extrapolate to WMAP frequencies, this aver- 
age flux translates into a flux smaller than 1 mJy in the W 
band. We thus reach a different conclusion about the relative 
importance of radio and infrared sources for White & Ma- 
jumdar (2004). This disagreemet is due to the fact that Lin 
et al. (2009) focus on galaxy clusters and also their extrap- 
olation is more accurate since their frequencies are closer to 
the W band. 



8.3 Contamination from lensed background 
sources 

An interesting mechanism capable of increasing the flux in- 
side a galaxy cluster with respect to the mean flux is lensing. 
Lensing of background galaxies by foreground clusters can 
create magnifications of a factor 10 or more. Big magnifi- 
cations occur far more often around galaxy clusters than in 
non-cluster regions. Also, in the mm to sub-mm band an in- 
teresting phenomenon occurs, the negative k-correction. A 
typical starburst galaxy will be seen with almost the same 
intensity at rcdshifts 2 and 5. This means that if one is able 
to reach a sensitivity enough to detect redshift 2 galaxies at 
these frequencies, the same sensitivity will allow one to see 
all these galaxies at even higher rcdshifts. Lensing of mm 
and sub-mm sources has been discussed in the past (a do- 
tailed analysis can be found in for instance Perrotta ct al. 
2003, see also Paciga et al. 2009) but from a statistical point 
of view and focusing on the number counts. 

In this section we focus on the possible contamination 
by lensed galaxies of the SZE. We use a toy model to quan- 
tify the level of contamination introduced by lensed galaxies. 
In our model we consider a cluster at a redshift of 0.15 (rep- 
resentative of the average redshift of our sample). Low red- 
shift clusters have larger lensing cross sections. This means 
that the lensing effect will be smaller for higher redshift clus- 
ters. The mass of the cluster is 8 x IO^^'Mq and is simulated 
by superposing three NFW profiles with different centers. 



The superposition of the NFW halos adds some substruc- 
ture to the model cluster but to first order it has spherical 
symmetry. In this case, lower masses will produce smaller 
magnifications. For this cluster we compute the magnificat 
tion map assuming a redshift of the source of z = 2. Ob- 
servations of mm and sub-mm sources with SCUBA at 850 
microns (Clements et al. 2008) and MAMBO at 1.2 mm 
(Bertoldi et al. 2007) show that the peak of the distribution 
is around 2 = 2. Given the redshift of our toy cluster, the 
magnification map depends weakly on the redshift of the 
sources. We use the rmmber counts derived with SCUBA 
(Coppin et al. 2006) and simulate a background of sources 
at redshift 2. The sources are distributed randomly and uni- 
formly. We make 100 different realizations and for each one 
we obtain their lensed counterpart. Finally, we compare the 
lensed with the unlensed images. We take discs of increasing 
radius centered on the cluster and compute the integrated 
fiux for each disc in the two images (lensed and unlensed). 
The result is shown in figure 13. The unlensed curves cor- 
respond to the dashed lines and they grow with the square 
of the disc radius (as expected). The lensed integrated fiux 
(dotted lines) behaves almost the same with the exception 
of a few realizations (« 7%) where significant magnifica- 
tion occurs. The maximum magnification occurs around the 
critical curves. Also, we found that a number of realizations 
(w 10%) suffer the contrary effect (demagnification of the 
integrated flux) caused by a large magnification of an under- 
dense (or empty) region. With this toy model we can con- 
clude that in the case of SCUBA sources, the lensing effect 
can increase significantly the flux in the cluster compared 
with the SZE flux in only a few percent of the cases. Also, 
by taking the average of a sample of clusters, magniflca- 
tion and demagnification effects compensate each other (at 
least partially) making the average lensing contribution even 
smaller. A study combining SCUBA (850 fim) and MAMBO 
(1200 /Km) sources revealed that the typical spectral index 
for these sources is a = 2 (Greve et al. 2008) where S cc . 
This implies that at lower frequencies the lensing effect is ex- 
pected to be even less relevant as there should be a smaller 
density of bright background sources (and hence, it would 
be more unlikely to have a source magnified). For example, 
from the number counts given in Bertoldi et al (2007), the 
density of sources at 1.2 mm is about an order of magni- 
tude less than the density of SCUBA sources at 1 mJy. For 
experiments operating at even lower frequencies (like the At- 
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Figure 13. Predicted contamination from lensing for an illustra- 
tive model. The curves show the integrated flux from background 
sources in a disc of radius K centered in the lensing cluster. For 
a uniform distribution of background sources the flux grows as 
i?^. Dashed lines show the integrated flux for the unlensed field 
of view. Dotted lines arc the corresponding flux for the Icnscd 
field of view. Each line corresponds to a different realization. The 
cluster is located at rcdshift 0.1-5 and it has a mass of 8x lO-'^Af,;:;. 
The background sources were assumed to be at redshift 2. A typ- 
ical SZE flux is of the order of a few hundred mjy at the SCUBA 
wavelength (850 /im). 

acama Cosmology Telescope, or ACT, Kosowsky 2006) the 
average impact of lensed sources should be even smaller with 
only a few percent of clusters showing a significant contam- 
ination to their SZE flux. The same conclusion applies to 
WMAP froquoncios where the contamination from lensing 
effect should be negligible. 



9 SZE PLUS POINT SOURCE MODEL 

In this section we try to fit the SZE profile considering a 
model that includes a point source (or PS) at the center of 
the cluster. Since the PS contributes as a positive AT and 
the SZE as a negative one, a partial cancellation of the later 
signal might occur. 

Lin et al. (2009) observed a sami)lc of 139 radio se- 
lected sources inside galaxy clusters. Tliesc sources were ob- 
served using the Very Large Array (or VLA) at frequencies 
between 4.9 and 43 GHz (note that the highest frequency 
matches the frequency band of the Q channel in WMAP). 
They found that most radio sources in clusters lay near the 
center of their host cluster and that a significant amount 
of the sources in clusters show a flat or even inverted spec- 
trum that would increase their flux at higher frequencies 
compared with the predicted flux using standard frequency 
extrapolations. We compute the average flux of their sample 
at the highest frequency (43 GHz) and find this to be 9.2 
mJy, or 12.4 mJy if we consider only the flatter spectrum 
unresolved and core sources. Some of the fluxes at 43 GHz 
given in Lin et al. (2009) were upper limits on the flux. Also, 
the minimum flux in their sample is a few mJy. Wc made a 
flt to their sample (see figure 14) using a simple power law 
and recomputed the average flux extrapolating down to 1 
mJy. The average flux for the model was 10.4 mJy. These 



numbers give us an idea that we should expect on average a 
PS contribution of about 10 mJy per cluster in the Q band. 
This valuable information can be used to set an additional 
constrain in our model. Also in Lin et al (2009) , the authors 
extrapolated the fluxes observed at lower frequencies to 90 
GHz. Two extrapolations were done using different assump- 
tions for the spectral index. The average flux of the sources 
in their table 3 at 90 GHz is 4.4 and 6.2 mJy for the two 
extrapolations. We should keep in mind however, that this 
extrapolation is highly uncertain and should be used with 
caution. Nevertheless, it is useful to establish that we can 
expect contributions of the order of a few mJy from radio 
sources at frequencies of 90 GHz (roughly, the W band in 
WMAP). 

Regarding infrared sources, the situation is less clear as 
there are very few observations of infrared sources in galax:y 
clusters. As noted earlier, from the results of Zemcov et al. 

(2007) using SCUBA data, one can infer a contribution of 
less than 1 mJy in galaxy clusters from thermal dust emis- 
sion. Other mechanisms (or the presence of colder gas in 
galaxies) could boost this flux to higher limits but this re- 
mains unclear so far. 

An interesting result can be found in Granett et al. 

(2008) . In this paper the authors stack CMB regions around 
50 supervoids and 50 superclustcrs found in the Sloan Dig- 
ital Sky Survey (or SDSS, Adclman-McCarthy et al. 2008). 
They found a colder average signal around the voids and a 
hotter average signal round the superclustcrs. Their inter- 
pretation of the results is that the colder and warmer signals 
are due to the integrated Sachs- Wolfe effect (Sachs & Wolfe 
1967). This is certainly one of the possible explanations but 
one could argue as well that the colder region is explained 
by a lack of point sources in the supervoid region and the 
warmer region by an excess of point sources in those regions. 
If the emission from point sources in these regions is domi- 
nated by mm sources with spectrum (Grevc et al. 2008), 
their frequency dependence would be virtually indistinguish- 
able from the ISW effect. Also, in Granett et al. (2008), the 
signal from superclusters should also be sensitive to the neg- 
ative contribution from the SZE. Instead, a positive signal 
is observed. Whether this positive signal is due to the ISW 
or point sources in the superclusters is an issue that needs 
to be resolved in the future for instance by observing the 
same regions at higher frequencies. 



10 DISCUSSION 

A model with a steeper profile (like the ADOS considered 
here) is capable of simultaneously reproducing the X-ray 
brightness and SZE profiles better /3-model. Lieu et al. 
(2006) predicted a large SZE signal in their sample of clus- 
ters. This is possibly due to the fact that they used an 
isothermal /3-model to fit the X-ray emission. We can argue 
that part of the discrepancies between the expected SZE 
signal considered in the past and the apparently smaller sig- 
nal seen by WMAP can be explained by considering models 
for the gas distribution that are steeper than the /3-model. 
This conclusion is in agreement with recent findings based 
on high resolution X-ray data in galax;y clusters that seem 
to suggest a steeper profile than previously thought (Neu- 
mann 2005, Vikhhnin et al. 2005, 2006, Croston et al. 2008). 
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Figure 14. Radio sources from Lin et al. (2009). The solid line 
shows a histogram of the radio sources observed in galaxy clusters 
by Lin et al (2009). The mean flux of these sources is 9.2 mjy in 
the Q band. The dotted line represents a model for the number 
counts which extrapolates the fluxes down to 1 mJy. The mean 
flux for this model is 10.4 iriJy. Lower minimum boundaries for 
the flux will result in a lower average flux. 

Atrio-Barandela et aL (2008) also found evidence of a steep- 
ening profile using WMAP data and a cluster sample similar 
to ours but using a different analysis. In their paper they 
concluded that the gas profile was more consistent with a 
NFW-like profile than a /3-model profile. Mroczkowski et 
al. (2009) has recently completed a study combining X-ray 
CHANDRA observations with SZE high resolution interfcro- 
metric measurements obtained with the Sunyaev-Zel'dovich 
Array (or SZA). In their analysis, they used a model with 
steeper profile than the 5-model profile for the pressure and 
they also found good agreement between the X-ray and SZE 
observations. The long baselines of the SZA also permit the 
detection of radio galaxies inside the clusters and thus an 
estimate of their contamination level. Mroczkowski et al. 
(2009) found that in their small sample of clusters, the ra^ 
dio source contamination was at the level of a few mJy at 
31 GHz. These kind of studies, were high resolution X-ray 
and SZE are combined in the same analysis, will be able to 
determine a more accurate model for the gas distribution. 
The interferometric data used in Mroczkowski et al. (2009) 
does not allow a study of the gas on scales larger than that 
set by the shortest baselines of the SZA. Future data from 
the Atacama Cosmology Telescope (ACT, Kosowski 2004) 
or the South pole Telescope (SPT, Ruhl et al. 2004) will 
complement the interferometric studies and extend the an- 
gular scale range. Also, the Planck satellite will provide a list 
of SZE measurements in many hundreds of clusters. Most of 
the clusters used in the present work (if not all) will be ob- 
served by Planck. Even though Planck will not be able to 
resolve most clusters, it will be useful to determine the to- 
tal flux of each cluster which can then be directly compared 
with the X-ray flux. This kind of comparison will allow a 
better determination of the gas contents in clusters. 

Perhaps one of the most revealing papers in the lit- 
erature about the incompatibility of the /3-modcl and X- 
ray-|-SZE data can be found in Hallman et al. (2007). In 
this paper the authors use simulated data of X-ray and SZE 
observations and conclude that "there is an inconsistency 



Figure 15. One dimensional profiles (solid lines) compared with 
three /3-models. Prom left to right, difference between models and 
X-ray RASS profile, Q, V and W WMAP profiles. The dashed 
line represents a /3 = 2/3 model similar to Mi but with p = 10. 
The dotted line is like the previous model but with /? = 0.5 and 
the dot-dashed line is for the same model but with /3 = 0.8 

between X-ray and SZE fitted model parameters that leads 
to a bias in deduced values of the Compton parameter and 
gas mass when using isothermal p-models - X-ray and SZE 
radial profiles from the same galaxy cluster can not be fit- 
ted by identical isothermal (3 -models" , Using our models we 
also find that it is hard to fit the X-ray and SZE profiles 
simultaneously using /3-models. Prom the previous discus- 
sion it seems reasonable to assume that higher values for 
f3 might agree better with the SZE measurements. We find 
that this is true but at the expense of failing to predict the 
X-ray measurements well. In figure 15 we show an example 
where we compare the average profiles measured by ROSAT 
and WMAP with models where we change the exponent 13. 
Higher values of f) seem to fit better the WMAP SZE pro- 
files (when the models are subject to the constraint that the 
total X-ray flux matches the measured X-ray flux). These 
models however, seem to under-predict significantly the X- 
ray emission in the wings of the proflle. A similar behavior 
is observed when other parameters in the /3-model, like the 
core radius or the parameter p are varied. Reducing the core 
radius (or p) brings the SZE down but at the expense of fail- 
ing to reproduce the wings in the X-ray profile. Although the 
/3-models do not seem to agree well with both observations 
(X-ray ad SZE) we find that in general, /3-models with large 
values for (3 perform better (than models with smaller (3) 
when compared to both data sets simultaneously. A model 
like the ADOS which is intrinsically steeper (or other similar 
ones) is able to reproduce even better the X-ray and SZE 
profiles simultaneously. 

When adopting our concordance model (M7 in table 
1), we still found evidence for missing SZE in WMAP data. 
We attempt to explain this with a model for the contami- 
nation from point sources (radio and infrared) inside galaxy 
clusters. We estimated that infra-red sources should give a 
very small contribution in flux of the order of 1 mJy or less 
at WMAP frequencies. Radio sources, on the other hand 
could contribute with several mJy in all WMAP bands al- 
though this estimation is generally derived from extrapo- 
lations made from lower frequency observations. In the Q 
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Figure 16. One dimensional profiles (solid lines) compared with 
the ADOS M7 model (dashed line). Prom left to right, X-ray RASS 
profile, Q, V and W WMAP profiles. The dotted Unes in the Q, V 
and W bands show the ADOS My concordance model minus the 
contribution of a point source at the center of the cluster with 
fluxes of 16, 26 and IS mjy in Q, V and W bands respectively. 

band a direct meEisurement of the average flux from point 
sources in a sample of clusters was made by Lin et al (2009). 
Radio sources arc expected to contribute about 10 mJy on 
average in this band. A fit combining the ADOS model and a 
point source contamination of 16 mJy seems to render a god 
fit to the data in the Q band (see figure 16). At higher fre- 
quencies it is still not clear whether the difference between 
the model and the data is caused by point sources, a differ- 
ent gas profile or other unidentified mechanism. Given the 
actual constraints on point source contamination (about 10 
mJy in the Q band and a several mJy in the W band) it 
may not be possible to reconcile the galaxy cluster X-ray 
and SZE observations without invoking more exotic or elab- 
orated models for the gas distribution in galaxy clusters. 
Some of our models were able to fit well the WMAP SZE 
profiles and simultaneously be consistent with the expected 
level of contamination from point sources in Q, V and W 
bands. These models though, seem to under-prcdict the X- 
ray emission from clusters at intermediate radii. The models 
that perform better at these radii seem to over-predict the 
SZE signal in WMAP. In this case, estimates of point source 
contamination provide some help to reconcile the model with 
the observations. In figure 16 we show one of those models 
where the point source contamination is 16, 26 and 18 mJy 
in Q, 1/ and W bands respectively. A /3- model with a good 
fit to the X.-ray profile would require point source contam- 
ination several times higher than these values which would 
contradict the measurements of Lin et al. (2009). Another in- 
teresting conclusion that can be derived from figure 16 is the 
fact that the point source contamination seems to be largest 
in the V band, in contradiction with previous expectations 
that extrapolate the radio emision to larger frequencies and 
the infrared emission to smaller ones. This might be an in- 
dication that the cluster model in figure 16 still needs to 
be improved but it also leaves the door open to future sur- 
prises in the frequency range 60-70 GHz. Data from Planck 
will help to answer this question in the incoming years. 

There axe other possibilities not explored in this paper 
to explain a possible lack of SZE signal. One is the local 



dumpiness of the gas around smaller dark matter sub-halos 
in the galaxy cluster. This dumpiness would increase the 
amount of X-rays produced with the same amount of gas. 
Hence, a fit to the X-ray surface brightness with a smooth 
model (like the /3-model or ADOS models) will again over- 
predict the amount of gas responsible for the X-ray emission 
and therefore the SZE produced by the same gas. A simi- 
lar over-prediction will occur when shocks and other non- 
thermal phenomena occur in galaxy clusters that boost the 
X-ray brightness surface. CHANDRA and XMM observa- 
tions reveal that some clusters show these features (Marke- 
vitch et al. 2000, McNamara et al. 2001, Fabian et al. 2003, 
Gioia ct al. 2004, Ma et al 2009). Lin et al (2004) discuss also 
another interesting scenario in which preheating of the gas 
in clusters occurs at high redshift [z ~ 2). This mechanism 
was studied earlier by Kaiser (1991), and Evrard & Henry 
(1991). Such heating reduces the SZE flux by a factor of 2-3 
on scales of a few arcminutes. 

One other interesting conclusion can be derived from 
our results. Earlier works have predicted and/or measured 
the existence of tight scaling relations between X-ray and 
SZE derived quantities (McCarthy ct al. 2003, Morandi et 
al. 2007, Bonamente et al. 2008, Yuan ct al. 2009). The 
use of scaling relations has proved to be useful in the past 
to derive other quantities such as the cluster temperature. 
They arc a powerful tool to investigate the physical prop- 
erties of the clusters and their evolution in redshift. Using 
our concordance model as the best guess for the underlying 
gas distribution, we can predict the total SZE flux from each 
cluster and compare it with the observed X-ray flux. To first 
order, we should expect a clear correlation between the two 
since they both depend on the square of the distance. The 
X-ray depends on the luminosity distance and the SZE on 
the angular diameter distance and for low redshifts both are 
very similar. In figure 17 we show the correlation for the 
M7 model. As expected, a clear correlation exists between 
the two quantities. The solid line represents the following 
scaling law; 



Flux = 45 X Sl^ (14) 



where Flux is the SZE fiux in mJy at 90 GHz and Sx is the 
X-ray flux in units of 10^^^ erg/s/cm^ and in the energy 
band (0.1-2.4) keV. The scatter around this scaling law is 
probably smaller than the real scatter in clusters but the 
scaling law itself can be considered representative as it was 
derived from the average properties of the cluster sample. 

Our results have profound implications for ongoing and 
future SZE surveys. Most of the predictions made in the 
past regarding the expected number of SZE detections for 
different experiments were made under the assumption that 
the gas profile can be well described by a /3-model. Our 
findings suggest that this assumption is not valid and that 
a steeper gas profile would be more adequate. A steeper 
profile subject to the constraint that it has to be compatible 
with X-ray observations results in a smaller amount of SZE. 
Future data like ACT, SPT or planck will answer this and 
other interesting questions. 
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Figure 17. Predicted SZE flux vs measured X-ray flux. The pre- 
dicted SZE was computed using the concordance model (Mr). 
The soUd Une represents a scaUng law (see text). 

11 CONCLUSIONS 

By stacking WMAP 5 year data wc arc able to sec the av- 
erage SZE signal from a large sample of galaxy clusters. 
These clusters axe X-ray selected. A comparison of models 
with the average X-ray profile and the average SZE profile 
allows us to narrow down the range of possible models. Stan- 
dard ,5-modcls arc harder to reconcile with both data sets. A 
steeper profile, like the one suggested by Ascasibar & Diego 
(2008) render a better fit to the data. None of the mod- 
els considered is able to perfectly fit both data sets. Some 
error can be allowed in the fits, in particular in the X-ray 
profile, where non-thermal mechanisms or dumpiness in the 
gas can boost the X-ray emissivity without requiring more 
gas, we can find a concordance model that predicts a sig- 
nal close to the observed ones (in X-rays and SZE) and also 
leaves room for point source contamination at roughly the 
level expected from radio and infrared extrapolations. We 
conclude that the X-ray and SZE measurements are consis- 
tent once a careful modeling of both data sets and a realistic 
level of point source contamination are taken into account. 
Based on the predictions made by the concordance model, 
we compare the SZE and X-ray fiuxes and find a tight corre- 
lation between them. These kinds of correlations have been 
predicted in the past but never observed with a large sample 
of clusters. Future SZE data like Planck, and also high reso- 
lution experiments like ACT or SPT, will help answer some 
of the questions raised in this paper and further develop this 
new and exciting field. 
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